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Background: The severe acute respiratory syndrome (SARS) virus has undergone mutations in its receptor-binding 


domain. 


Results: We used biochemical, functional, and crystallographic methods to investigate these mutations. 


Conclusion: These mutations were viral adaptations to either the human or palm civet receptor. 
Significance: This research elucidates detailed mechanisms of host receptor adaptation by the SARS virus and can help predict 
and monitor future evolution of the SARS virus in animals. 


The severe acute respiratory syndrome coronavirus (SARS- 
CoV) from palm civets has twice evolved the capacity to infect 
humans by gaining binding affinity for human receptor angio- 
tensin-converting enzyme 2 (ACE2). Numerous mutations have 
been identified in the receptor-binding domain (RBD) of differ- 
ent SARS-CoV strains isolated from humans or civets. Why 
these mutations were naturally selected or how SARS-CoV 
evolved to adapt to different host receptors has been poorly 
understood, presenting evolutionary and epidemic conun- 
drums. In this study, we investigated the impact of these muta- 
tions on receptor recognition, an important determinant of 
SARS-CoV infection and pathogenesis. Using a combination of 
biochemical, functional, and crystallographic approaches, we 
elucidated the molecular and structural mechanisms of each of 
these naturally selected RBD mutations. These mutations either 
strengthen favorable interactions or reduce unfavorable inter- 
actions with two virus-binding hot spots on ACE2, and by doing 
so, they enhance viral interactions with either human (hACE2) 
or civet (CACE2) ACE2. Therefore, these mutations were viral 
adaptations to either hACE2 or cACE2. To corroborate the 
above analysis, we designed and characterized two optimized 
RBDs. The human-optimized RBD contains all of the hACE2- 
adapted residues (Phe-442, Phe-472, Asn-479, Asp-480, and 
Thr-487) and possesses exceptionally high affinity for hACE2 
but relative low affinity for cACE2. The civet-optimized RBD 
contains all of the cACE2-adapted residues (Tyr-442, Pro-472, 
Arg-479, Gly-480, and Thr-487) and possesses exceptionally 
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high affinity for cACE2 and also substantial affinity for hACE2. 
These results not only illustrate the detailed mechanisms of host 
receptor adaptation by SARS-CoV but also provide a molecular 
and structural basis for tracking future SARS-CoV evolution in 
animals. 


Many viruses, particularly RNA viruses, can rapidly adapt to 
changes in their hosts due to the high mutation rates of these 
viruses. Such viral adaptation to hosts imposes huge threats to 
human health because it allows viruses to emerge from animals 
to infect humans. Understanding viral adaptation to hosts 
would allow humans to track viral evolution that may lead to 
viral epidemics, but current knowledge of the mechanisms of 
viral adaptation to hosts is limited. This study investigates the 
mechanisms of host receptor adaptation by the severe acute 
respiratory syndrome coronavirus (SARS-CoV), the agent of 
the SARS epidemic, and provides a molecular and structural 
basis for monitoring future SARS-CoV evolution in animals. 

SARS broke out as a worldwide epidemic in 2002-2003, 
causing >8000 human infections with a 10% fatality rate (1-4). 
SARS briefly recurred in 2003-2004, causing sporadic human 
infections with mild symptoms and no human-to-human 
transmission (5-7). Since the SARS epidemic, various SARS- 
CoV strains have been isolated from humans and palm civets 
(see Fig. 1A). The prototypic human strain, hTor02, was iso- 
lated during the 2002-2003 SARS epidemic (8). A strain com- 
mon to humans and civets, hcGd03, was isolated during the 
2003-2004 sporadic SARS infections (6, 7). Strain cSz02 was 
isolated from civets in animal markets in southern China dur- 
ing the SARS outbreak (5). Strain cCHb05 was isolated from wild 
civets in central China several years after the SARS outbreak 
(9). In addition, human strain hHae08 was isolated after adap- 
tation of SARS-CoV to human airway epithelial cells in vitro 
(10). Because of the extremely high genetic similarity (~99%) 
between human and civet SARS-CoV, it is believed that civets 


3 The abbreviations used are: SARS-CoV, severe acute respiratory syndrome 
coronavirus; SLCV, SARS-like coronavirus; ACE2, angiotensin-converting 
enzyme 2; hACE2, human ACE2; cACE2, civet ACE2; RBD, receptor-binding 
domain; RBM, receptor-binding motif. 
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played a direct role in transmitting SARS-CoV to humans and 
were responsible for both of the human SARS infections (5, 11, 
12). Although SARS-like coronaviruses (SLCVs) have been iso- 
lated from bats, the genetic difference between SLCVs and 
human or civet SARS-CoV is much greater than that between 
human and civet SARS-CoV (13, 14). Moreover, SLCVs cannot 
infect human or civet cells, and SARS-CoV cannot infect bats 
(15). Because of the direct role of civets in the past SARS epi- 
demic, this study focused on SARS-CoV strains isolated from 
humans and civets. 

Receptor recognition is the first step in viral infection of 
host cells and one of the most important determinants of 
viral infectivity and pathogenesis. Although many other host 
and viral factors can also affect the efficiency of infection and 
replication of a virus in a specific host, these factors come 
into play only after a virus has bound to a receptor on the cell 
membrane. A trimeric spike protein on the SARS-CoV enve- 
lope mediates viral entry into host cells. It first binds to its 
host receptor, angiotensin-converting enzyme 2 (ACE2) 
(16), and subsequently fuses host and viral membranes (17— 
19). A defined receptor-binding domain (RBD) on the SARS- 
CoV spike protein is sufficient for high-affinity binding to 
ACE2 (20-22). The RBD/ACE2 binding affinity has been 
established as a key element of SARS-CoV pathogenesis and 
cross-species infections (10, 11, 23-25). The crystal struc- 
ture of hTor02 RBD complexed with human ACE2 (hACE2) 
showed that the RBD contains two subdomains, a core struc- 
ture and a receptor-binding motif (RBM) that directly binds 
ACE2 (see Fig. 1B) (26). Using hTor02 RBD as a reference, 
numerous naturally selected mutations have been identified 
in the RBDs of different SARS-CoV strains isolated from 
humans or civets. The majority of these mutations fall in the 
RBM region (see Fig. 1, A and B). To date, research has 
focused on two of these mutations, N479K and T4878 (7, 11, 
24, 27, 28). Biochemical studies showed that both mutations 
N479K and T4875 significantly decrease the RBD/hACE2 
binding affinity (24, 28). Structural studies revealed that 
mutation N479K introduces an energetically unfavorable 
positive charge at the RBD/hACE2 interface, whereas muta- 
tion T487S removes an energetically favorable hydrophobic 
interaction with hACE2 Lys-353 (see Fig. 1C) (26, 27). The 
molecular and structural mechanisms of all of the other 
RBM mutations are poorly understood. It is not known why 
any of these RBM mutations were naturally selected or which 
RBM mutations were viral adaptations to which host. 

Here, we investigated the role of each of the RBM muta- 
tions in viral interactions with hACE2 and civet ACE2 
(cACE2) using a combination of biochemical, functional and 
crystallographic methods. Our results show that all of the 
above RBM mutations were viral adaptations to two “virus- 
binding hot spots” on ACE2 and suggest that these muta- 
tions were naturally selected because they enhance viral 
interactions with either hACE2 or cACE2. These detailed 
mechanisms allow us to understand the events that led to the 
past SARS epidemic and to monitor the events that may lead 
to future SARS epidemics. 
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EXPERIMENTAL PROCEDURES 


Protein Purification, Crystallization, and Structure Deter- 
mination—SARS-CoV RBDs (from different viral strains) and 
ACE2 peptidase domains (human, civet, or chimeric) were 
expressed and purified as described previously (26, 27). In brief, 
the proteins containing an N-terminal honeybee melittin signal 
peptide and a C-terminal His tag were expressed in Sf9 insect 
cells using the Bac-to-Bac system (Invitrogen). The proteins 
were harvested from Sf9 cell supernatants and purified sequen- 
tially on a nickel-nitrilotriacetic acid column and a Superdex 
200 gel filtration column (GE Healthcare). To purify the RBD- 
ACE2 complex, ACE2 was incubated with excess RBD, and 
then the complex was purified by gel filtration chromatogra- 
phy. RBD-ACE2 crystals were grown in sitting drops at room 
temperature over wells containing 100 mm Tris (pH 8.5), 
20-24% PEG 6000, and 100 mm NaCl. Crystals were soaked 
briefly in 100 mm Tris (pH 8.5), 30% PEG 6000, 100 mm NaCl, 
and 30% ethylene glycol before being flash-frozen in liquid 
nitrogen. X-ray diffraction data were collected at Argonne Pho- 
ton Source beamlines 19-ID and 24-ID and Advanced Light 
Source beamline 4.2.2. Each structure was determined by 
molecular replacement using the structure of hACE2 com- 
plexed with hTor02 RBD as the search model (Protein Data 
Bank code 2AJF). Data and refinement statistics are shown in 
Table 1. The software used for data processing and structure 
refinement is also listed in Table 1. 

RBD/hACE2 Binding Assays—Surface plasmon resonance 
assays using a Biacore 2000 system were carried out as 
described previously (29, 30). In brief, to measure the binding 
affinities between mutant viral RBDs and wild-type ACE2, 
ACE2 was immobilized on a C5 sensor chip through direct 
covalent coupling via amine groups, and soluble RBDs were 
introduced at a flow rate of 20 l/min at different concentra- 
tions (62.5, 125, and 250 nm). To measure the binding affinities 
between mutant ACE2 and wild-type viral RBDs, RBDs were 
immobilized on the sensor chip, and ACE2 was flowed through. 

Pseudotyped Viral Infection Assays—Infection was assayed 
using murine leukemia viruses expressing B-galactosidase and 
pseudotyped with SARS-CoV spike protein as described previ- 
ously (29). In brief, to prepare pseudotyped viruses, HEK293T 
cells were cotransfected with spike protein-encoding pcDNA3.1, 
p3240 expressing murine leukemia virus gag and pol genes, and 
a murine leukemia virus B-galactosidase-transducing vector 
(pBAG) (31). Viral supernatants were harvested and inoculated 
on HEK293T cells expressing ACE2 in pcDNA3.1. Infection 
efficiency was quantified by measuring B-galactosidase 
activity. The spike proteins and ACE2 contained a C-termi- 
nal C9 tag and HA tag, respectively. Hence, the concentra- 
tions of the spike protein packaged in pseudotyped viruses 
and of ACE2 expressed on the HEK293T cell surface were 
detected by Western blotting using anti-C9 and anti-HA 
antibodies, respectively, and quantified using Image) soft- 
ware (version 1.6) (supplemental Fig. $1). Compared with 
the infection experiments using fluorescence assays (32), the 
B-galactosidase assays used in this study are not as sensitive 
but are more stable and repeatable. 
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FIGURE 1. Interface between SARS-CoV RBD and hACE2. A, list of mutations in the RBMs of various SARS-CoV strains. Five representative existent strains and 
two predicted future strains are defined in the Introduction. B, overall structure of the hTor02 RBD-hACE2 complex (Protein Data Bank code 2AJF). hACE2 is in 
green, and hTor02 RBD is in cyan (core) and red (RBM). RBM residues that underwent mutations are displayed. C, detailed structure of the hTor02 RBD/hACE2 
interface. hACE2 residues are in green, SARS-CoV residues that underwent mutations are in magenta, and SARS-CoV residues that played significant roles in the 


mutations are in cyan. 


RESULTS AND DISCUSSIONS 


Two Virus-binding Hot Spots on hACE2—All of the RBM 
mutations cluster around two virus-binding hot spots on 
hACE2, hot spot-31 and hot spot-353, which center on Lys-31 
and Lys-353, respectively (Fig. 1C). Hot spot-353 is recognized 
by both SARS-CoV and human NL63 coronavirus, despite the 
markedly different structures of their RBDs (29, 30). Hot 
spot-31 is recognized by SARS-CoV only but not by NL63 coro- 
navirus. The two hot spots have both structural similarities and 
dissimilarities. In unbound hACE2, both Lys-31 and Lys-353 
point away from the protein and into solution (33). Upon 
SARS-CoV binding, both residues fold back, become buried in 
hydrophobic environments, and form salt bridges with hACE2 
Glu-35 and Asp-38, respectively (26, 27). Because of the hydro- 
phobic environments, formation of the salt bridges provides 
significant amounts of energy to the virus/receptor interac- 
tions. The difference between the two hot spots is that, upon 
salt bridge formation, Lys-353 is in a relatively extended con- 
formation that is energetically favorable, whereas Lys-31 is in a 
strained conformation that is energetically unfavorable. In this 
study, we examined the structure and function of hot spot-31 
by mutagenesis using surface plasmon resonance Biacore 
assays and pseudotyped viral infection assays (Fig. 2, A and B). 
Both assays showed that mutation K31A did not have a signif- 
icant impact on RBD/hACE2 binding (p > 0.10 (¢ test) for both 
K,, and infection). In contrast, mutation K353A significantly 
decreased RBD/hACE2 binding affinity (v < 0.01 for both K, 
and infection). These results suggest that there was a significant 
net energy gain from salt bridge formation at hot spot-353. In 
contrast, there was no significant net energy gain from salt 
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bridge formation at hot spot-31, as the energy gain from salt 
bridge formation was cancelled out by loss from the strained 
conformation of Lys-31. 

To further understand the structure and function of hot 
spot-31, we investigated mutation K31T in hACE2 based on the 
observation that cACE2, containing Thr-31, is an efficient 
receptor for hTor02 (Fig. 2C). Both Biacore and pseudotyped 
viral infection assays showed that mutation K31T significantly 
increased RBD/hACE2 binding affinity (py < 0.001 for both K, 
and infection) (Fig. 2, A and B). To investigate the structural 
basis for mutation K31T, we determined the crystal structure of 
a human-civet chimeric ACE2 complexed with hTor02 RBD 
(Table 1). The chimeric ACE2 contained SARS-CoV-binding 
residues from cACE2 (Fig. 2C) and other residues from hACE2. 
We previously established this chimeric ACE2 as a useful alter- 
native to cACE2 in crystallographic analysis of SARS-CoV/ 
cACE2 interactions by showing that the chimeric ACE2 had the 
same receptor activities as cACE2 but the same crystallo- 
graphic activities as hACE2 (27). The structure of the hTor02 
RBD-chimeric ACE2 complex showed that a hydrogen bond 
between chimeric ACE2 Thr-31 and RBD Tyr-442 strengthens 
a hydrophobic network around hot spot-31 at the virus/recep- 
tor interface (Fig. 2D). Therefore, despite its insignificant net 
energy contribution, hot spot-31 plays an important role in 
RBD/ACE2 interactions. Consequently, upon introduction to a 
new species, SARS-CoV RBD may need to adapt to both hot 
spot-31 and hot spot-353 on ACE2. 

Molecular Mechanisms of RBM Mutations in SARS-CoV— 
To investigate the molecular mechanisms of the RBM muta- 
tions, we introduced seven single mutations to hTor02 RBD, 
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FIGURE 2. Structures and functions of two virus-binding hot spots on hACE2. A, surface plasmon resonance Biacore analysis of the binding interactions 
between hTor02 RBD and wild-type or mutant hACE2. hTor02 RBD was immobilized, and wild-type or mutant ACE2 was flowed through. Each experiment was 
repeated six times at three different protein concentrations. The corresponding S.E. values are shown. B, pseudotyped viral infection assays of the interactions 
between hTor02 spike protein and wild-type or mutant hACE2. Retroviral murine leukemia viruses expressing B-galactosidase and pseudotyped with hTor02 
spike protein were used to infect HEK293T cells expressing wild-type or mutant hACE2. Infection efficiency of pseudotyped viruses was measured by B-galac- 
tosidase assays and normalized against the infection efficiency in cells expressing wild-type hACE2. Each experiment was repeated six times. The correspond- 
ing S.E. values are shown. C, use of a human-civet chimeric ACE2 in crystallographic studies of RBD/cACE2 interactions. The chimeric ACE2 contains SARS-CoV- 
binding residues from cACE2 and other residues from hACE2. The chimeric ACE2 has the same receptor activities as cCACE2 but the same crystallographic 
activities as hACE2 (27). SARS-CoV-binding residues that differ between hACE2 and cACE2 are shown. D, structure of the interface between hTor02 RBD and the 


hTor02 RBD 


chimeric ACE2. 


TABLE 1 
Crystallographic data collection and refinement statistics 


X-ray diffraction data were processed using HKL2000 (35). Programs CNS (36) and CCP4 REFMAC (37) were used for structure refinement. Structural illustrations were 


made using POVscript (38). r.m.s.d., root mean square deviation. 


Complex 


hTor02 RBD-chimeric ACE2 


hOptimize RBD-hACE2 


cOptimize RBD-chimeric ACE2 cOptimize RBD-hACE2 


Data 
Space group P2, P2 
Cell constants a = 81.6, b = 119.5, 
c= 113.4A; B = 92.5° 


1 
a = 814, b = 118.3, 
c= 111.9 A; B = 93.15° 


P2, P2, 
a= 81.0, b = 120.0, a = 81.2, b = 119.3, 
c= 112.2 A; B = 93.15° c= 113.2 A; B = 92.25° 


Resolution (A) 50-3.0 50-2.9 50-3.0 50-3.0 

Roy (last shell) (%) 9.8 (61.6) 7.5 (33.4) 7.0 (51.6) 9.2 (48.0) 

Observed reflections 157,444 223,385 122,747 296,564 

Unique reflections 42,738 47,509 42,870 44,284 

Completeness (last shell) (%) 97.7 (95.1) 93.4 (59.1) 91.7 (59.9) 95.1 (67.0) 

I/o (last shell) 6.3 (1.2) 11.4 (3.2) 16.5 (1.6) 20.2 (2.2) 
Refinement 

Reon (Rive) (%) 23.9 (29.2) 22.6 (28.3) 23.9 (28.5) 23.0 (27.8) 

Ree reflections (%) 5 5 5 5 

Bond length r.m.s.d. (A) 0.012 0.010 0.011 0.010 

Bond angle r.m.s.d. 1.361° 1.252° 1.329° 1.257° 

Chiral r.m.s.d. (A?) 0.092 0.090 0.086 0.083 


Y442F, L472P, L472F, N479K, N479R, D480G, and T4875, cor- 
responding to all of the RBM mutations identified in human 
and civet SARS-CoV strains (Fig. 1A). We expressed and puri- 
fied each of these mutant RBDs and measured their binding 
affinities for hACE2 and cACE2 using Biacore assays (Fig. 3, A 
and B). The results showed that hTor02 RBD bound to hACE2 
and cACE2 with similar K, values (py > 0.50). The affinity of 
hTor02 RBD for hACE2 was increased by Y442F and L472F 
(p < 0.01); unchanged by N479R (p > 0.10); and decreased by 
L472P, D480G, N479K, and T4875 (p < 0.02) (Fig. 3A). The 
affinity of hTor02 RBD for cACE2 was increased by N479R and 
D480G (p < 0.01); unchanged by L472P and N479K (p > 0.50); 
and decreased by L472F, Y442F, and T487S (p < 0.01) (Fig. 3B). 
In sum, among the RBM residues that have been identified in 
different SARS-CoV strains, Phe-442, Phe-472, Asn-479/Arg- 
479, Asp-480, and Thr-487 enhanced RBD/hACE2 interac- 
tions, whereas Tyr-442, Pro-472/Leu-472, Arg-479, Gly-480 
and Thr-487 enhanced RBD/cACE2 interactions. 
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The above results suggest that combinations of the RBM resi- 
dues that enhance RBD/hACE2 or RBD/cACE2 interactions can 
potentially produce new RBDs with exceptionally high affinity for 
hACE2 or cACE2, respectively. To test this hypothesis, we 
designed and constructed two such new RBDs: hOptimize RBD 
containing Phe-442, Phe-472, Asn-479, Asp-480, and Thr-487 and 
cOptimize RBD containing Tyr-442, Pro-472, Arg-479, Gly-480, 
and Thr-487 (Fig. 1A). We expressed and purified each of the 
RBDs. Biacore assays confirmed that compared with hTor02 RBD, 
hOptimize and cOptimize RBDs had higher affinities for hACE2 
and cACE2, respectively (p < 0.01) (Fig. 3, A and B). On the other 
hand, whereas hOptimize RBD had low affinity for cACE2, cOp- 
timize RBD had substantial affinity for hACE2. 

Structural Mechanisms of RBM Mutations in SARS-CoV— 
To investigate the structural mechanisms of the RBM muta- 
tions, we determined the following three crystal structures: 
hOptimize RBD complexed with hACE2, cOptimize RBD com- 
plexed with the chimeric ACE2, and hOptimize RBD complexed 
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FIGURE 3. Molecular mechanisms of RBM mutations in SARS-CoV. A, surface plasmon resonance Biacore analysis of the binding interactions between hACE2 
and wild-type or mutant hTor02 RBD. hACE2 was immobilized, and wild-type or mutant hTor02 RBD was flowed through. B, surface plasmon resonance Biacore 
analysis of the binding interactions between cACE2 and wild-type or mutant hTor02 RBD. cACE2 was immobilized, and wild-type or mutant hTor02 RBD was 
flowed through. C, pseudotyped viral infection assays of the interactions between different spike proteins and hACE2. D, pseudotyped viral infection assays of 


the interactions between different spike proteins and cACE2. 
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FIGURE 4. Structural mechanisms of RBM mutations in SARS-CoV. A, structure of the interface between hACE2 and hOptimize RBD. B, structure of the 
interface between the chimeric ACE2 and cOptimize RBD. C, structure of the interface between hACE2 and cOptimize RBD. D, superimposed structures of the 
interfaces between the chimeric ACE2 and cOptimize RBD (colored) and between hACE2 and hTor02 (white). 


with the chimeric ACE2 (Fig. 4 and Table 1). Together, these struc- 
tures show that naturally selected RBM mutations either strength- 
ened favorable interactions or reduced unfavorable interactions 
with the two virus-binding hot spots on ACE2. The structural 
details of each of the RBM mutations are as follow. 

RBD residue 442 interacts with hot spot-31 on ACE2. At the 
RBD/hACE2 interfaces, the hydroxyl group of Tyr-442 has par- 
tial steric clash with the alkyl side chain of hACE2 Lys-31 (Fig. 
1C). Mutation Y442F removes this hydroxyl group, partially 
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releases the structural strain of Lys-31 (Fig. 4A), and increases 
the RBD/hACE2 binding affinity. As discussed above, at the 
RBD/chimeric ACE2 interfaces, the hydroxyl group of Tyr-442 
forms a critical hydrogen bond with chimeric ACE2 Thr-31 
(Figs. 2D and 4B). Mutation Y442F removes this hydrogen 
bond, weakens the hydrophobic network around hot spot-31, 
and decreases the RBD/cACE2 binding affinity. Therefore, 
Phe-442 and Tyr-442 were viral adaptations to hACE2 and 
cACE2, respectively (Fig. 5). 
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FIGURE 5. Summary of host receptor adaptation by SARS-CoV. Listed are 
adaptations of RBM residues to hACE2 or cACE2. Arrows point from less well 
adapted residues to better adapted residues. Double arrows connect equally 
well adapted residues. 


RBD residue 472 interacts with hot spot-31 on ACE2. At the 
RBD/hACE2 interfaces, Leu-472 forms favorable hydrophobic 
interactions with hACE2 Met-82 and Leu-79 and strengthens 
the hydrophobic network around hot spot-31 (Fig. 1C). Muta- 
tion L472F further strengthens but mutation L472P weakens 
RBD/hACE2 interactions (Fig. 4, A and C). At the RBD/chime- 
ric ACE2 interfaces, Phe-472 would have partial steric clash 
with the hydroxyl group of chimeric ACE2 Thr-82 (Figs. 2D and 
4B), and hence, mutation L472F decreases the RBD/cACE2 
binding affinity. Additionally, neither Leu-472 nor Pro-472 has 
significant contact with Thr-82 (Figs. 2D and 4B), and hence, 
mutation L472P has no significant effects on the RBD/cACE2 
interactions. Therefore, Phe-472 and Pro-472/Leu-472 were 
viral adaptations to hACE2 and cACE2, respectively (Fig. 5). 

RBD residue 479 interacts with both hot spots on ACE2. As 
shown previously, at the RBD/hACE2 interfaces, mutation 
N479K introduces a positive charge that can compete with 
hACE2 Lys-31 for hACE2 Glu-35 as a salt bridge partner (26, 
27), and hence, it decreases the RBD/hACE2 binding affinity. 
Unexpectedly, unlike N479K, mutation N479R has no signifi- 
cant effects on the RBD/hACE2 binding affinity even though it 
does introduce a positive charge. This is because Arg-479 and 
hACE2 Lys-353 can both form a salt bridge with hACE2 Asp-38 
(Fig. 4C). At the RBD/chimeric ACE2 interfaces, Thr-31 cannot 
form a salt bridge with chimeric ACE2 Glu-35 (Figs. 2D and 
4B), and hence, Glu-35 is available to form a salt bridge with 
Lys-479. On the other hand, Lys-479 has high conformational 
entropy and stacks poorly with the hydrophobic network 
around hot spot-31 (27). Consequently, mutation N479K has 
no significant effects on the RBD/cACE2 binding affinity. In 
contrast, Arg-479 can form a strong bifurcated salt bridge 
with Glu-35 and is also a much better hydrophobic stacker 
than Lys-479 (Fig. 4B). Thus, mutation N479R increases the 
RBD/cACE2 binding affinity. Therefore, Asn-479/Arg-479 
and Arg-479 were viral adaptations to hACE2 and cACE2, 
respectively (Fig. 5). 

RBD residue 480 interacts with hot spot-353 on ACE2. At the 
RBD/hACE2 interfaces, RBD Tyr-436 forms a hydrogen bond 
with hACE2 Asp-38 and strengthens the Lys-353—Asp-38 salt 
bridge (Figs. 1C and 4D). The alkyl side chain of Asp-480 forms 
a hydrophobic interaction with Tyr-436 and supports the Tyr- 
436—Asp-38 hydrogen bond. Consequently, mutation D480G 
weakens the above interactions and decreases the RBD/hACE2 
binding affinity. At the RBD/chimeric ACE2 interfaces, chime- 
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ric ACE2 Glu-38, which has a longer side chain than hACE2 
Asp-38, moves back to keep the salt bridge with Lys-353. 
Accordingly, Tyr-436 also needs to move to keep the hydrogen 
bond with Glu-38 (Fig. 4D). Mutation D480G facilitates the 
movement of Tyr-436 and increases the RBD/cACE2 binding 
affinity. Therefore, Asp-480 and Gly-480 were viral adaptations 
to hACE2 and cACE2, respectively (Fig. 5). 

RBD residue 487 interacts with hot spot-353 on hACE2. As 
shown previously, at the RBD/hACE2 interfaces, the methyl 
group of Thr-487 supports the Lys-353—Asp-38 salt bridge (Fig. 
1C) (26, 27), and hence, mutation T487S decreases the RBD/ 
hACE2 binding affinity. At the RBD/chimeric ACE2 interfaces, 
the methyl group of Thr-487 also supports the Lys-353—Glu-38 
salt bridge (Fig. 4B), and thus, mutation T4875 decreases the 
RBD/cACE2 binding affinity. Therefore, Thr-487 was a viral 
adaptation to both hACE2 and cACE2 (Fig. 5). 

Evaluation of Host Cell Infections Mediated by Different 
SARS-CoV Spike Proteins—To further investigate the interac- 
tions between different RBDs and ACE2 from humans or civets, 
we measured the spike protein-mediated infection of hACE2- or 
cACE2-expressing cells using pseudotyped viral infection assays 
(Fig. 3, Cand D). Compared with Biacore assays that directly mea- 
sure binding interactions between monomeric RBD and mono- 
meric ACE2 in solution, pseudotyped infection assays measure 
binding interactions between trimeric spike proteins on viral sur- 
faces and monomeric ACE2 on target cell surfaces. In addition, 
pseudotyped infection assays require not only receptor binding 
but also membrane fusion and post-fusion steps for reporter gene 
expression. Hence, pseudotyped infection assays reflect more 
truthfully than Biacore assays on how the RBM mutations may 
impact SARS-CoV infectivity and pathogenesis. 

The spike proteins used in the assays share the same 
sequences as the hTor02 spike protein except that their RBDs 
contain different combinations of RBM residues (Fig. 1A). 
hOptimize spike protein mediated a higher level of infection 
of hACE2-expressing cells than all other spike proteins (p < 
0.01) because its residues at positions 442, 472, 479, 480, and 
487 were all best adapted to hACE2 (Fig. 3C). Particularly, 
hOptimize spike protein-pseudotyped viruses infected 
hACE2-expressing cells about twice as efficiently as hTor02 
spike protein-pseudotyped viruses and ~4-fold as efficiently 
as cSz02 spike protein-pseudotyped viruses. Additionally, 
hOptimize spike protein mediated a lower level of infection 
of cACE2-expressing cells than all other spike proteins 
except hHae08 (p < 0.01) because its residues at four of the 
above positions, except position 487, were poorly adapted to 
cACE2. On the other hand, cOptimize spike protein medi- 
ated a higher level of infection of cACE2-expressing cells 
than all other spike proteins (yp < 0.01) because its residues at 
the above five positions were all best adapted to cACE2 (Fig. 
3D). Particularly, cOptimize spike protein-pseudotyped 
viruses infected cACE2-expressing cells about twice as effi- 
ciently as hTor02 spike protein-pseudotyped viruses and 
~7-fold as efficiently as hHae08 spike protein-pseudotyped 
viruses. cOptimize spike protein also mediated a substantial 
level of infection of hACE2-expressing cells, lower than 
hTor02 but significantly higher than hcGd03 (p < 0.01), 
because its residues at positions 479 and 487 were well 
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adapted to hACE2, whereas its residues at the other three 
positions were poorly adapted to hACE2. Overall, the results 
from the pseudotyped infection assays correlate well with 
those from Biacore data and structural analysis. 

Monitoring Future SARS-CoV Evolution in Animals—SARS- 
CoV evolution in animals needs to be carefully monitored as a 
preventive measure for potential future SARS epidemics. SARS- 
CoV was found in wild civets in 2005 (9) and can potentially cross 
the species barrier to infect humans again, as it did as least twice 
before. Novel coronaviruses have been identified in bats (13, 14), 
and some of these bat coronaviruses may be transmitted to 
humans either directly or through some intermediate animal hosts 
(11). SARS-CoV may also be misused as a bioweapon to infect 
humans. Regardless of the circumstances under which SARS-CoV 
re-emerges, understanding the mechanisms of its host receptor 
adaptation will be crucial for dealing with its re-emergence. 

In this study, we have identified RBM residues that are best 
adapted to hACE2 or cACE2, and hence, these residues can be 
predictors of future SARS-CoV evolution in animals. The elu- 
cidated receptor recognition mechanisms by SARS-CoV can 
allow us to evaluate the potential epidemical consequences of a 
future SARS-CoV strain simply by examining its genomic 
sequence. For example, if a future SARS-CoV strain is isolated 
from civets or humans, we will be able to evaluate the binding 
affinity between its RBD and cACE2 or hACE2 based on its 
RBM sequences. More specifically, if a future SARS-CoV strain 
contains a majority or all of the cACE2-adapted RBM residues, 
it can potentially have high infectivity for civet cells; if it con- 
tains a majority or all of the hACE2-adapted RBM residues, it 
can potentially have high infectivity for human cells. Moreover, 
the detailed interactions between SARS-CoV RBD and the two 
virus-binding hot spots on ACE2 as elucidated in this study will 
help us understand and evaluate new RBM mutations in future 
SARS-CoV strains that are different from the previously iden- 
tified RBM mutations. The results from this study can also con- 
tribute to the monitoring and control of SLCVs in bats by pro- 
viding insights into how these bat viruses may interact with the 
ACE2 molecules from bats, civets, humans, or other animals 
(27, 34). It should be noted that our study did not include rep- 
lication-competent viruses or animal models, and thus, future 
studies will be required to establish the correlation between 
optimized RBDs and infectivities of live SARS-CoV strains. 
To conclude, given that SARS-CoV may still be evolving in 
civets and that SLCVs are harbored by bats, it remains pos- 
sible for SARS-CoV to re-emerge into the human population 
through civets or other animals, but we can be better pre- 
pared to prevent or handle future outbreaks of the virus with 
improved understanding of the mechanisms for its host 
receptor adaptation. 
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